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The radiative capture cross sections for the 3 He(af,y) 7 Be and 3 H(af,y) 7 Li reactions are calculated in the fully 
microscopic fermionic molecular dynamics approach using a realistic effective interaction that reproduces the 
nucleon-nucleon scattering data. At large distances bound and scattering states are described by antisym- 
metrized products of 4 He and 3 He/ 3 H ground states. At short distances the many-body Hilbert space is extended 
with additional many-body wave functions needed to represent polarized clusters and shell-model-like config- 
urations. Properties of the bound states are described well, as are the scattering phase shifts. The calculated 
S factor for the 3 He(of,y) 7 Be reaction agrees very well with recent experimental data in both absolute normal- 
ization and energy dependence. In the case of the 3 H(a,y) 7 Li reaction the calculated S -factor is larger than 
available experimental data by about 15%. 
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The 3 He(a,y) 7 Be reaction is one of the key reactions in 
the solar proton-proton chains 111 0] . It competes with the 
3 He( 3 He,2/?) 4 He reaction and therefore determines the pro- 
duction of 7 Be and 8 B neutrinos in the ppll and ppIII branches. 
For a long time the experimental situation regarding the cap- 
ture cross section was not clear due to conflicting experimen- 
tal results [1]. In recent years the capture cross section has 
been remeasured at the Weizmann institute OH, by the LUNA 
Collaboration an, by the Seattle group |6|], and by the 
ERNA Collaboration |7|] now providing consistent high pre- 
cision data. Nevertheless, it is still not possible to reach the 
low energies relevant in solar burning and the data have to be 
extrapolated with the help of models. A careful analysis of 
the new data sets and a discussion of the extrapolation and its 
uncertainties is given in Ref . 12] . 

The first attempts to model the capture cross sections were 
done by using an external capture model li, @] where only 
the asymptotic form of the bound and scattering state wave 
functions enters, neglecting the behavior of the wave function 
at short distances. In potential models like, e.g., Ref. [10] 
the wave functions are described by two pointlike clusters in- 
teracting via an effective nucleus-nucleus potential which is 
adjusted to give reasonable properties for the bound states 
and the scattering phase shifts. In the framework of the mi- 
croscopic cluster model, e.g., Refs. flll-[l41. the system is 
described by antisymmetrized wave functions of two clus- 
ters. One has to solve for the relative motion of the clus- 
ters by using resonating group or generator coordinate meth- 
ods. In these microscopic models phenomenological nucleon- 
nucleon interactions are used. Like in the potential models 
these interactions are tuned to reproduce certain properties of 
bound and scattering states within the restricted cluster model 
space. There have been attempts lEHIH] to go beyond the 
single-channel approximation by including the 6 Li+/? chan- 
nel, but such enlarged model spaces require again modifica- 
tions of the phenomenological interaction. 

Predictive power is expected from ab initio methods which 
use realistic interactions that reproduce the nucleon-nucleon 



scattering data and the deuteron properties. Solving the many- 
body problem with realistic interactions is hard, as very large 
model spaces are required and up to now consistent ab ini- 
tio reaction calculations have been possible only for single 
nucleon projectiles Ml 1711 . The 3 He(a,y) 7 Be reaction was 
studied in hybrid approaches, where asymptotic normaliza- 
tion coefficients calculated from 7 Be bound state wave func- 
tions using variational Monte Carlo 111 811 and the no-core shell- 
model lll9ll were combined with conventional potential mod- 
els. None of these calculations is successful in describing both 
the normalization and the energy dependence of the capture 
cross section data. 

In this Letter, we present the first ab initio type calculation 
of the 3 He(a,y) 7 Be and 3 H(a,y) 7 Li capture cross sections. We 
describe consistently bound and scattering states starting from 
a realistic effective interaction derived in the unitary correla- 
tion operator method. The fermionic molecular dynamics ap- 
proach is used to create many-body wave functions that cap- 
ture the relevant physics in the interaction region. Frozen clus- 
ter configurations with 4 He and 3 He/ 3 H ground states are used 
at large distances. 

The effective interaction is derived from the realistic Ar- 
gonne VI 8 interaction J201, that reproduces the deuteron 
properties and the nucleon-nucleon scattering phase shifts. 
The interaction is transformed into a phase- shift equivalent 

by using 
El HI 

tral and tensor correlations are incorporated explicitly. In this 
work we use UCOM correlation functions that are derived 
from a Hamiltonian evolved using the similarity renormaliza- 
tion group (SRG) as described in Ref. J22I1 with a flow param- 
eter a = 0.20 fm 4 , corresponding to a soft cut-off A= 1.5 fm" 1 . 
As shown in Ref. 1E2I1 . no-core shell-model calculations using 
the two-body UCOM(SRG) interaction are able to reproduce 
the binding energies of triton, 4 He and 7 Li. 

Fermionic Molecular Dynamics (FMD) is a microscopic 
many-body approach that has been used successfully for nu- 
clear structure studies of nuclei in the p and sd shell. See 



low-momentum interaction by using the unitary correlation 
operator method (UCOM) JUl |22|] where short-range cen- 
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FIG. 1. (color online). Cuts through the density distributions of in- 
trinsic basis states. Left: Selected polarized configurations obtained 
in variation after angular momentum and parity projection for 3/2", 
7/2~ and 1/2" states. Right: Frozen configurations where only the 
cluster distance is varied. 



I23U25I1 for some recent applications and ll26ll for a general 
discussion. FMD is based on intrinsic many-body basis states 
that are Slater determinants 
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The complex parameters bk encode the mean positions and 
momenta of the wave packets. The width parameters are 
variational and can be different for each nucleon. The spin 
can assume any direction, and isospin & is ±1/2. The wave 
packet basis is very flexible and contains harmonic oscillator 
shell-model and Brink- type cluster states as special cases. 

To restore the symmetries of the Hamiltonian the intrinsic 
basis states \Q) are projected on parity, angular momentum 
and total linear momentum 
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so that the wave function factorizes into the internal part and 
the center-of-mass motion given by a plane wave. 

In general the intrinsic states have no axial symmetry and 
K is not a good quantum number. Linear dependent combina- 
tions among the different ^-projections have to be removed. 
This is done numerically and introduces a small ambiguity in 
the size of the model space. We will exploit this ambiguity 
later to fine-tune the 7 Be and 7 Li binding energies. 

All bound and scattering states are represented by using a 
set of intrinsic states |<2 W ) 
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Proper boundary conditions for bound and scattering states are 
imposed by using the microscopic /^-matrix approach devel- 
oped by the Brussels group 11271 12811 . 

At large distances the 7 Be and 7 Li wave functions consist 
of 4 He and 3 He/ 3 H clusters in their ground states interacting 
via the Coulomb interaction only. The relative motion of these 
frozen clusters is therefore given by Whittaker and Coulomb 
functions for bound and scattering states, respectively. Micro- 
scopically we describe these cluster configurations with FMD 



Slater determinants where the clusters are put at a distance R. 
The wave functions of the individual clusters are obtained by 
variation in the FMD model space. 

In the interaction region the nuclear interaction will polar- 
ize the clusters. To include these polarization effects we ex- 
tend the model space with additional FMD basis states ob- 
tained by variation after projection (VAP) on spin-parity 3/2", 
1/2", 7/2", and 5/2" as well as on l/2 + , 3/2 + , and 5/2 + . 
The square radius of the intrinsic state is used as a constraint 
to generate configurations corresponding to cluster distances 
from 1 to 5 fm. Together with the frozen configurations 
that extend to distances slightly beyond the channel radius 
a =12 fm, we have about 50 intrinsic basis states to represent 
the inner part of the wave function. Density distributions of 
typical frozen and polarized basis states are shown in Fig.[T] 

When the model space is restricted to frozen configurations 
the 3/2" and 1/2" states in 7 Be are bound by only 240 and 
10 keV respectively. The FMD VAP configurations are there- 
fore essential to a get a good description of the bound states. 
As mentioned the numerical elimination of linear dependent 
states in the ^-mixing procedure introduces a small ambiguity 
in the model space size that translates into an ambiguity in the 
binding energy of about 150 keV. As the reaction cross sec- 
tion depends very sensitively on phase space we exploit this 
ambiguity to tune the centroid of the 3/2" and 1/2" bound 
state energies to the experimental value. The calculated split- 
ting between the bound states is too small compared to the 
data. However, the total cross section essentially depends only 
on the centroid energy, whereas the branching ratio slightly 
changes with the splitting. The bound state properties for 
7 Be and 7 Li are summarized in Table J] The charge radii and 
quadrupole moments test the tail of the wave functions and 
agree reasonably well with experiment. 

In Fig.0 we show the phase shifts for scattering in the S- 
and Z)-wave channels. As for the bound states, the addition 
of polarized configurations to the model space significantly 
changes the results and leads to a good agreement with the 
available data JHH. 

The capture cross section for the 3 He(of,y) 7 Be reaction 
is calculated by using the many-body scattering and bound 
eigenstates of the Hamiltonian. In the energy range up to 
2.5 MeV, it has been shown lfl8ll that only dipole transitions 



TABLE I. Calculated and experimental bound state properties. Ener- 
gies with respect to the 4 He- 3 He and 4 He- 3 H thresholds, respectively. 
Experimental charge radii are from Refs. J29|,[3(3], the 7 Li quadrupole 
moment from Ref . Bill . 
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FIG. 2. (Color online). 4 He- 3 He scattering phase shifts. Dashed 
lines show results using only frozen configurations, solid lines are 
obtained with the full FMD model space. The calculated D-wave 
phase shifts lie on top of each other. Experimental results are from 
Refs. H3,[H. 



from the S - and Z)-wave scattering states have to be consid- 
ered. The obtained S factor is shown in Fig. [3] together with 
the experimental data. Our results are in good agreement with 
the recent measurements regarding both the absolute normal- 
ization and the energy dependence. The extrapolated zero- 
energy S factor is S 34 (0) = 0.593 keVb. 

As our model successfully describes the 3 He(of,y) 7 Be re- 
action, it should also do well for the isospin mirror reaction 
3 H(a,y) 7 Li. As shown in Fig. 0] we observe a good agreement 
for the energy dependence of the S -factor but find that the 
absolute normalization is about 15% larger than the data by 
Brune, Kavanagh, and Rolfs 113411 , 

In summary our calculations are able to describe consis- 
tently the bound state properties, and the scattering phase 
shifts as well as the normalization and energy dependence of 
the 3 He(a,y) 7 Be capture cross section. Our results deviate 



from the correlation between the ground state quadrupole mo- 
ment and zero-energy S factor found in cluster models using 
phenomenological interactions (3, H] . Our approach differs 
in two main aspects from those earlier studies. First, we use 
a well defined effective interaction that describes the nucleon- 
nucleon scattering data. In contrast to phenomenological ef- 
fective interactions the UCOM interaction has a pronounced 
momentum dependence and a longer range due to the explic- 
itly included pion exchange, a feature that turns out to be im- 
portant for the low energy scattering solutions. Second, our 
model space is larger than in the cluster model. Additional 
FMD basis states in the interaction region describe polarized 
clusters and shell-model-like configurations. Although they 
are only a small admixture in the full wave functions they are 
essential to describe the bound state properties as well as the 
scattering phase shifts. 

The results can also be studied in terms of overlap func- 
tions that are obtained by mapping the microscopic many- 
body wave functions onto the relative wave function of two 
pointlike nuclei in the resonating group formalism. In Fig. \5\ 
we show the overlap functions for the l/2 + scattering state 
at E cm = 50keV and the 3/2" bound state. The nodes in the 
overlap functions reflect the antisymmetrization between the 
clusters. We also show the dipole strength calculated with 
these overlap functions. It reproduces the dipole matrix ele- 
ment calculated with the microscopic wave functions within 
2%. Comparing with the dipole strength obtained from the 
Coulomb and Whittaker functions matched at the channel ra- 
dius we observe sizable differences up to distances of about 
9 fm. This indicates that the assumption of predominant ex- 
ternal capture at low energies is not that well satisfied. 

Future calculations should investigate the role of three- 
body forces. It is expected that low-momentum three body 
forces would increase the splitting between the 3/2" and 1 /2" 
states but would have a minor effect on the centroid energy. 
Furthermore more detailed wave functions could be used. In 
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FIG. 3. (Color online). The astrophysical S factor for the 
3 He(of,y) 7 Be reaction. The FMD result is given by the solid line. 
Recent experimental data 0-0] are shown as dark colored symbols 
and older data (it] as light symbols. 



FIG. 4. (Color online). The astrophysical S factor for the 3 H(af,y) 7 Li 
reaction. The FMD result is given by the solid line. Most recent 
experimental data are shown as dark symbols and older data as light 
symbols (II34I1 and references therein). 
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FIG. 5. (Color online). On the left: 7 Be overlap functions for a low 
energy l/2 + scattering state and the 3/2~ bound state (solid lines). 
Coulomb and Whittaker functions matched at the channel radius 
(dashed lines). On the right: Dipole strength calculated with over- 
lap functions (solid line) and with Coulomb and Whittaker functions 
(dashed line). 



the FMD approach it is difficult to describe long-range tensor 
correlations explicitly, so that the absolute binding energies 
are underestimated, although the binding energy with respect 
to the cluster threshold is in very good agreement with no- 
core shell-model results. Nevertheless, we expect that such 
improvements will not change the capture cross sections sig- 
nificantly as important properties like phase shifts of the scat- 
tering states, binding energy with respect to the cluster thresh- 
old, asymptotic behavior of the bound state wave functions as 
tested by charge radius and quadrupole moment, and proper 
treatment of antisymmetrization are already well described in 
the present calculation. 
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